PUPT-2297 



Lectures on Holographic Superfluidity and Superconductivity 

C. P. Herzog 

Department of Physics, Princeton University, Princeton, NJ 08544, USA 

(Dated: March 26, 2009) 

Four lectures on holography and the AdS/CFT correspondence applied to condensed matter 
systems.^ The first lecture introduces the concept of a quantum phase transition. The second 
lecture discusses linear response theory and Ward identities. The third lecture presents 
transport coefficients derived from AdS/CFT that should be applicable in the quantum 
critical region associated to a quantum phase transition. The fourth lecture builds in the 
physics of a superconducting or superfluid phase transition to the simple holographic model 
of the third lecture. 

I. INTRODUCTION: QUANTUM PHASE TRANSITIONS 

Spurred by the concrete proposal of refs. [1-3] for an AdS/CFT correspondence, there are some 
good reasons why holographic ideas have become so important in high energy theoretical physics 
over the last ten years. The first and perhaps most fundamental reason is that the AdS/CFT 
conjecture provides a definition of quantum gravity in a particular curved background space-time. 
The second is that AdS/CFT provides a tool for studying strongly interacting field theories. These 
lectures concern themselves with the second reason, but I will spend a paragraph on the first. 

Given the lack of alternative definitions of quantum gravity, the AdS/CFT conjecture is diffi- 
cult to prove, but the correspondence does give a definition of type IIB string theory in a fixed 
ten dimensional background and by extension of type IIB supergravity. Recall that the original 
conjecture posits an equivalence between type IIB string theory in the space-time ^^5.5 x and 
the maximally supersymmetric (SUSY) SU(A'^) Yang-Mills theory in 3-M dimensions. (In our nota- 
tion, AdS^ is five dimensional anti-de Sitter space and is a five dimensional sphere.) Yang-Mills 
theory, at least in principle, can be simulated on a computer as the continuum limit of a lattice 
theory. The low energy limit of type IIB string theory is type IIB supergravity and the correspon- 
dence must also yield a quantum theory of gravity. This line of reasoning has led to an improved 
understanding of black hole physics, including a tentative resolution of the black hole information 
paradox and a better understanding of black hole entropy. 

Nearly as fundamental and no less exciting is the prospect of using AdS/CFT to understand 
strongly interacting field theories by mapping them to classical gravity. To see how the correspon- 
dence can be used as a tool, recall that the interaction strength of maximally SUSY Yang-Mills 
theory is described by the 't Hooft coupling A = gy^^N. Through the AdS/CFT correspondence 
A = {L/lg)'^ where L is the radius of curvature of AdS^ (and the S^), and Ig is a length scale that 
sets the tension of the type IIB strings.^ Strings are also characterized by a coupling constant. 



^ These lectures were given at the Spring Scliool on Superstring Tlieory and Related Topics, 23 - 31 March, 2009, 

ICTP, Miramare, Trieste, Italy. 
^ Tlie tension is conventionally defined as l/(27r£s)- 
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gs, that describes their UkeHhood to break. The AdS/CFT dictionary relates the string coupHng 

constant to the gauge theory coupHng via Angyyi = gs- In the double scaling limit where N —>■ oo 
while A is kept large and fixed, string theory is well approximated by classical gravity. Keeping 
is/L small means string theory is well approximated by gravity, while keeping gs small eliminates 
quantum effects. Using the AdS/CFT correspondence, enormous progress has been made over the 
last ten years in understanding the large A'^, large A limit of maximally SUSY Yang-Mills theory 
in 3+1 dimensions and its cousins. 

Before moving on to an extensive discussion of condensed matter systems, let me briefly mention 
that one of the most interesting ideas surrounding this flurry of activity mapping out the properties 
of strongly interacting SUSY field theories with gravitational duals is that we might learn something 
about quantum chromodynamics (QCD). At low energy scales, QCD is a quintessential example 
of a strongly interacting field theory. Consider temperatures slightly above the deconfinement 
transition, of the order of 200 MeV, where the baryons and mesons dissolve into a soup of strongly 
coupled quarks and gluons. Such a non-Abelian soup is probably not so qualitatively different 
from maximally SUSY Yang-Mills in the double scaling limit. Here is one example of an intriguing 
application of holographic techniques. Experiments at the relativistic heavy ion collider (RHIC) 
combined with hydrodynamic simulations suggest that the viscosity of the quark-gluon plasma is 
very low (see for example ref. [4]). In contrast, perturbative QCD techniques yield a large viscosity 
[5, 6], and lattice gauge theory requires a very difficult analytic continuation from the Euclidean 
theory to extract such a transport coefficient (see for example ref. [7]). AdS/CFT yields, for 
maximally SUSY Yang-Mills (and indeed for all its cousins in this double scaling limit) the low 
value rj/s = h/AirkB for the viscosity to entropy density ratio [8]. 

Putting aside QCD, in these lectures I will describe progress in applying holographic techniques 
to condensed matter systems. QCD is not the only useful strongly interacting field theory. Field 
theory has for a long time been a standard tool in a condensed matter theorist's toolbox. For 
example, near phase transitions, coherence lengths become long enough to allow a continuum 
description of a crystal lattice or otherwise discretized system of atoms and molecules. 

There are some structural and probabilistic reasons why QCD may not be the best candidate 
for an application of holographic techniques. Asymptotic freedom of QCD means at high energy 
scales the gravity dual will necessarily become increasingly stringy, and the correspondence loses 
some of its simplicity and power. While finding a holographic dual for QCD may be like finding a 
needle in the haystack of generalized AdS/CFT correspondences, the odds of finding a gravity dual 
to a condensed matter system appear, at least superficially, to be better. There are hundreds of 
thousands of pre-existing materials to consider. Moreover, using nano-lithography, optical lattices, 
and other experimental techniques, we may be able to engineer a material with a gravity dual. This 
last possibility raises the tantalizing prospect of better understanding quantum gravity through 
material science or atomic physics. 

An outline for the rest of these lectures is as follows: 

• Using the notion of a quantum phase transition, in the rest of the first lecture I will frame 
the connection between condensed matter systems and holography in a useful and hopeful 
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way. 

• The second lecture is a discussion of old and doubtless well known results in field theory. I 
have devoted a whole lecture to these results for a few reasons. The first is that it is much 
easier to understand what extra information AdS/CFT is giving us if we first understand 
the limitations of field theory. The second is that while the first, third and fourth lectures 
may not stand the test of time, the contents of this second lecture are true and probably 
very useful in other contexts. 

• In the third lecture, I will holographically compute field theory transport coefficients using 
a very simple gravitational action consisting of an Einstein-Hilbert and Maxwell term: 

S=^J d^x^iR -2A)-^J d^x^F^^F^'- . (1) 

I will connect these holographic results, in a qualitative way, to measurements of transport 
coefficients in graphene and high temperature superconductors. 

• In the last lecture, justifying the title of this lecture series, I will modify the gravitational 
action by adding an order parameter that will produce a superconducting or superfiuid 
phase transition. I will focus on the case where the order parameter is a scalar, but one 
could introduce a vector order parameter as well by promoting the Abelian F^i, in the action 
above to an SU(2) gauge field. 

A. Quantum Pheise Treinsitions 

The notion of a quantum phase transition in condensed matter systems provides our motivation 
for using AdS/CFT. A quantum phase transition is a phase transition between different phases of 
matter at T = 0. Such transitions can only be accessed by varying a physical parameter, such as a 
magnetic field or pressure, at T = 0. They are driven by quantum fluctuations associated with the 
Heisenberg uncertainty principle rather than by thermal fluctuations. We will be concerned with 
second order quantum phase transitions in this lecture. Much of the discussion here is drawn from 
ref. [9]. 

At r = but away from a quantum critical point, a system typically has an energy scale A 
perhaps associated with the energy difference between the ground and first excited state. Another 

important quantity is a coherence length ^ characterizing the length scale over which correlations 
in the system are lost. Let g be the physical parameter driving the quantum phase transition. At 
the quantum critical point we expect A to vanish and ^ to diverge, but not necessarily in the 
same way: 

A~(5-5cr, (2) 

i-ia-gcT". (3) 

The quantity z, relating the behavior A ~ is usually called the dynamical scaling exponent. At 
the quantum critical point, the system becomes invariant under the rescaling of time and distance. 
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t X^t and X — > Ax. Different z occur in different condensed matter systems. For example, 
z = 1 is common for spin systems, and we will see an example of such a system shortly. The 
case z = 1 is special because the quantum critical system typically has a Lorentz symmetry and 
the scaling becomes a part of a larger conformal symmetry group SO{d + 1,2) for a system in d 
spatial dimensions. These lectures will focus mostly on the z = 1 case because it is here that the 
AdS/CFT dictionary is most powerful and well developed. Another common and familiar value 
is z = 2. The free Schrodinger equation is invariant under z = 2 scalings, but there are other 
examples as well, e.g. Lifshitz theories. Generic, non-integer z are possible. 

Figure 1 shows a prototypical phase diagram for a system that undergoes a quantum phase 
transition. Here the physical parameter is a coupling g, and the quantum phase transition occurs 

9 = 9c and T = 0. At low temperatures, we imagine the system is in one of two phases 
well characterized by some order parameter (s). The solid blue lines in the phase diagram could 
be classical thermal phase transitions or softer cross-overs, depending on the dimensionality and 
nature of the system. The region between the dashed black lines is the quantum critical region 
(QCR). 
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FIG. 1: A typical phase diagram involving a second order quantum critical point. 

The usefulness of the notion of a quantum phase transition lies in a wished for ability to 
understand the system in the QCR. The QCR is characterized by the requirement that T be large 
compared to the dimensionally appropriate power of {g — gc)- It seems reasonable to expect that the 
effective scale invariant field theory valid at the critical point, now generalized to nonzero T, can be 
used to predict the behavior of the system in the QCR. (We can generalize this discussion, replacing 
T with some other external parameter or set of parameters — chemical potential, magnetic field, 
etc.) 
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B. The Qucintum Rotor 

The quantum rotor is a simple theoretical model that exhibits a quantum phase transition. The 
model is described by the Hamiltonian: 

H = 9jJ2Lf-jJ2ni-nj , (4) 

where we are summing over a lattice (of arbitrary dimension) indexed by i and where (ij) indicates 
a pair of nearest neighbor sites. Let hi be an A'^ component vector such that hf = 1. The operator 
Li is an angular momentum, and Lf is thus the kinetic energy term for this vector hi which lives 
on an — 1 dimensional sphere. Taking J > 0, the interaction term in H will prefer to align the 
hj. The kinetic energy, in contrast, is minimized by randomizing the hi such that (L?) = 0. (For 
more details about this model, see ref. [9].) 

The quantum rotor exhibits a quantum phase transition as we tune the value oi g. In the limit 
5 S> 1, the sites on the lattice decouple from one another, and the system can be solved exactly. 
In the ground state, the kinetic energy is minimized by taking (hi) = such that (L|) = 0. 
Correlations between different lattice sites die off exponentially with distance, 

(0|ni-n,|0)~e-l"*-^^l/« , (5) 

where ^ is the correlation length. The lowest energy excitation is a particle where a single lattice 
site has a nonzero (Lf), and this particle hops from site to site. There is an energy gap A+ ~ gj 
associated with this particle. Because an external field will tend to align the hi, the ground state 
in this limit is a paramagnet. 

In contrast, in the opposite limit g <^ 1, the system becomes magnetically ordered. It is 
energetically favorable that {hi) / and for all of the vectors to align: 

lini {0\hi-hj\0) = Ni . (6) 

\Xi—Xj\— too 

This alignment spontaneously breaks rotational symmetry, and there must be an associated mass- 
less Nambu-Goldstone boson. These massless bosons are spin waves, i.e. slow rotations in the 
direction of {hi). Although there is no energy gap associated with this continuum of excited spin 
wave states, one can define an energy scale A_ associated with the kinetic term of the spin waves 
— a spin stiffness. (High energy theorists might prefer the term pion decay constant.) 

It seems reasonable to infer that there is a quantum phase transition between these two different 
types of order for a critical value g = gc of the coupling. Referring to Figure 1, phase one for this 
model would be magnetically ordered while phase two is the quantum paramagnet. As we approach 
the critical point, it should be energetically easier for the spin waves to rotate more quickly in the 
magnetically ordered phase or for particle excitations to form in the paramagnetic phase. Thus, as 
g approaches gc, we expect the energy scale to vanish as a power of (g — gc). Also, as we move 
out of the paramagnetic phase, the correlation length ^ should diverge as the vectors hi align. 

The quantum rotor is more than a toy. On the experimental side, for two spatial dimensions and 
N = 3, ref. [9] argues that the Hamiltonian models two sheets of La2Cu04, the parent compound of 
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a high Tc superconductor I will discuss at greater length in a moment. On the theoretical side, the 
continuum limit of this model should be very familiar to field theorists. It is the 0(iV) nonlinear 
sigma model. The continuum Lagrangian takes the form 



subject to the constraint = 1. It is thus a model with dynamical exponent 2; = 1, as was 

promised for spin systems. Note that for condensed matter applications, the effective speed of 
light would typically be much less than the actual speed of light, c^s <C c. This continuum limit 
should become a better and better approximation close to the quantum critical point where the 
correlation length diverges and we can coarse grain the hi degrees of freedom. 
More loosely, we could soften the constraint and replace it with a scalar potential 



to control the size of the fluctuations, yielding the O(A^) vector model. At this point, I refer the 
reader to a standard field theory textbook such as ref. [10] for a more thorough treatment than I 
can provide here. In the renormalization group language, the 0(iV) vector model is known to flow 
to a strongly interacting Wilson-Fisher fixed point, which in our language is nothing other than a 
quantum critical point. 



Quantum phase transitions are believed to be important in describing superconducting- insulator 
transitions in thin metallic films, as is demonstrated pictorially by rotating Figure 2 ninety degrees 
counter-clockwise. The rotated diagram is meant to resemble closely Figure 1 where phase one is 
an insulator, phase two is a superconductor, and g corresponds to the thickness of the film. The 
insulating transition is a cross-over, while the superconducting transition might be of Kosterlitz- 
Thouless type. There exists a critical thickness for which the system reaches the quantum critical 
point at T = 0. 

One of the most exciting (and also controversial) prospects for the experimental relevance of 
quantum phase transitions is high temperature superconductivity. Consider the parent compound 
La2Cu04 of one of the classic high Tc superconductors, La2-xSr^Cu04. La2Cu04 is actually not 
a superconductor at all but an anti-ferromagnetic insulator at low temperatures. The physics of 
this layered compound is essentially two dimensional. The copper atoms are arranged in a square 
lattice on separated sheets with effectively one electron per unit cell. The spins of the electrons 
pair up in an anti-ferromagnetic order. 

To turn La2Cu04 into a superconductor, the compound can be doped with strontium which has 
the effect of removing one electron for every lanthanum atom replaced with strontium. Figure 3 is 
a phase diagram for La2-a;Sra;Cu04. Once the doping x becomes sufficiently large, the compound 
superconducts at low temperature. Introducing some vocabulary, the doping which yields the 
highest Tc is called the optimal doping; for this material, Xo ~ 3/20 yielding a Tc 40 K. When 




(7) 



V{n) = a\n{x)\'^ + /3\n{x)\'^ , 



(8) 



C. Quantum Critical Points in the Real World 
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FIG. 2: Resistivity of thin films of bismuth versus temperature. The different curves correspond to different 
thicknesses, varying from a 4.36 A film that becomes insulating at low temperatures, to a thicker 74.27 A 
film that becomes superconducting. The figure is reproduced from ref. [11]. 

X > Xo, the compound is referred to as over doped, while when x < Xq, the compound is called 
under doped. 

Over doped high Tc superconductors are better understood than their under doped counterparts. 
For temperatures T > Tc, the material behaves like a Fermi liquid where quasiparticle, electron-like 
degrees of freedom are effectively weakly coupled. Moreover, the phase transition seems to follow 
the BCS paradigm where the electrons form Cooper pairs as we lower the temperature below Tc. 
In contrast, in the under doped region, the effective degrees of freedom are believed to be strongly 
interacting. In one paradigm, the superconducting to normal phase transition involves disordering 
the phase of the condensate rather than breaking Cooper pairs, if it indeed makes sense to talk 
about quasiparticles at all in this regime. Because the electrons may remain in bound states in the 
normal, under doped region of the phase diagram, this region is sometimes called the pseudogap. 
For more details about these issues, the reader might try ref. [12]. 

Speculations about the relevance of a quantum phase transition are related to the dip in Tc at 
a doping of Xc = 1/8 and a possible connection between this dip in Tc and experimental evidence 
for so-called striped phases where spin and charge density waves break translational invariance at 
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FIG. 3: A cartoon phase diagram for a superconductor such as La2_a;Sr2.Cu04. AF stands for anti- 
ferromagnetic and SC for superconducting. 

distance scales of order of a few times the lattice spacing [13]. The conjecture is that we can add 
a third axis to our phase diagram corresponding to an extra control parameter g in some model 
Lagrangian for the system, as pictured in Figure 4. In this figure, the chemical potential n plays 
the role of doping. In the third direction, we may find a quantum critical point where the dip in 
Tc becomes more pronounced and reaches the T = plane. One might hope to gain theoretical 
control over the pseudogap region using the effective field theory of the quantum critical point. 

The reader may ask why we have invented this extra quantum critical point when, by our 
definition, there appear to be three perfectly good quantum critical points already present in the 
phase diagram of Figure 3, one of which is close to the under doped region we would like to 
understand. The answer is ultimately unsatisfactory and is indicative of the speculative nature of 
the last few paragraphs. Not all quantum phase transitions are created equal, and this putative 
fourth quantum critical point promises to be a little simpler and cleaner. Although there is little 
data available, of the original three critical points, the transition from the normal phase to the 
superconducting phase in the under doped region is probably first and not second order. Moreover it 
is believed to be disorder driven and thus involves breaking translation invariance at large distance 
scales.^ 



Thus far, I have tried to argue that quantum phase transitions are important in understanding 
superconducting insulator transitions in thin films and may be important in the physics of high Tc 



D. On the Role of AdS/CFT 



I would like to thank Subir Sachdev and Markus Mueller for discussion on this point. 
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FIG. 4: A third conjectural axis has been added to our phase diagram for a high Tc superconductor. This 
figure was taken from ref. [13]. 

superconductors. I also exhibited the quantum rotor, a model theoretical system which undergoes 
such a phase transition. 

In general, it can be difficult to describe a system at a strongly interacting quantum critical 
point. Weakly coupled effective degrees of freedom may be difficult to identify or not exist, as I 
sketched for the case of high Tc superconductors. The reader may ask, "Can't we always discretize 
the system and simulate it on a computer?" To avoid problems with oscillatory numerical integrals, 
lattice models are almost always formulated for computers in Euclidean time. For questions about 
equilibrium physics (with no chemical potential), the answer is often, "Yes, the lattice is good 
enough." However, if we want to ask questions about physics at nonzero density, about real 
time physics, about transport coefficients and response to perturbations, numerical lattice models 
require tricky and usually untrustworthy analytic continuations. 

Returning to the role of holography in this story, AdS / CFT provides a tool to study a class of 
strongly interacting field theories with Lorentz symmetry in d space-time dimensions by mapping 
them to classical gravity in d+l space-time dimensions. The correspondence is a very useful way of 
working out the equation of state, real time correlation functions and transport properties such as 
diffusion constants, conductivities, and viscosities. The ambitious program is to find an example of 
an AdS/CFT correspondence that describes a real world material. Less ambitiously, we may learn 
universal or semi-universal properties about a class of strongly interacting field theories. These are 
often field theories and questions for which AdS/CFT is our only calculational tool. 

More specifically, when I talk about an AdS/CFT correspondence that describes a real world 
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material, I mean given a material that undergoes a quantum phase transition, the gravity dual 
should provide a good effective description of the critical point and the critical region in Figure 1. 
We might even be able to model the thermal phase transitions away from the critical region, but 
eventually the microscopic degrees of freedom in the material will become important, and a field 
theory description will be less valuable. 

The reader may object that there could well be structural reasons why this program is doomed to 
fail. AdS/CFT correspondences typically involve some underlying supersymmetry while condensed 
matter systems do not. In mitigation, I note that introducing chemical potential and temperature 
breaks supersymmetry, that if we stay away from the T = and fx = limits, the physics may not 
be so different whether the underlying theory is supersymmetric or not. 

The reader may also complain that the restriction to z = 1 appears to be limiting given the many 
different types of scaling that appear in condensed matter. In response, I note there has been recent 
progress in extending AdS/CFT to z 1. For example refs. [14, 15] have conjectured a gravity dual 
for a theory with Schrodinger symmetry, i.e. the symmetry group of the free Schodinger equation. 
This group has z = 2. There is serious experimental motivation to understand strongly interacting 
systems with this symmetry group [16]. Consider a dilute gas of lithium-6 or potassium-40 atoms 
in an optical trap. The interaction strength between these fermionic atoms can be tuned with an 
external magnetic field. At a Feshbach resonance, the scattering length becomes larger than the 
system size, and these so-called fermions at unitarity obey an approximate Schrodinger symmetry. 
Although I will not address this question here, it is very interesting to ask whether AdS/CFT can 
say anything useful about these strongly interacting atomic systems. 

Not all condensed matter systems with z = 2 have the Schrodinger symmetry. Another possi- 
bility is a Lifshitz scaling symmetry, i.e. the symmetry group of a Lagrangian of the form 

£ = {dt<pf - K{V^f . (9) 

Ref. [17] presents a proposal for a gravity dual for a strongly interacting field theory with such a 
symmetry. Despite these recent advances, we shall focus henceforth on the z = 1 case with Lorentz 
symmetry. It is this case for which the AdS/CFT dictionary is most detailed and reliable. 
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II. FIELD THEORY FOR STRONGLY INTERACTING SYSTEMS 

In this lecture, I will set up a field theory framework to describe the response of a system 
at equilibrium to small perturbations. The framework allows us to relate two-point correlation 
functions to thermal and charge conductivities. I then will demonstrate how gauge and Lorentz 
invariance severely constrain the form of these two-point functions, and hence of the conductivities, 
through Ward identities.^ 

Consider the response of a system to the presence of weak external fields {0i(a;)} coupled to a 
set of operators {0^{x)}. The Hamiltonian H is modified by a term of the form 

dH = - J d'^x (t>i{t, x)d\t, x) . (10) 

A classic result from time dependent perturbation theory in quantum mechanics, which I leave 
as a worthwhile exercise for the reader, is that these external fields will produce a change in the 
expectation value of the operators of the form 

d{&{x)) = J d'^+^x'G%{x,x')^j{x') + 0{(t?) , (11) 

where 

Gl{x,x')=ie(t-t'){[d\x),&{x')]) (12) 

is the retarded Green's function. In systems with translation invariance, it is sensible to decompose 
the external potential into Fourier components. The Fourier transformed linear response result 
takes the simple form 

50\k) = Gl{k)~^^{k) + 0{<t>^) , (13) 
where the Fourier transform of the retarded Green's function is 

G'^{k) = j d'^+^xe-'''''G%{x,0) . (14) 

A. The Relation Between Green's Functions and Treinsport Coefficients 

To make this formalism more intuitive and physical, I will rephrase Ohm's Law in this language. 
Recall the statement of Ohm's Law for an electric field that is constant in space but oscillating in 
time with frequency to; the spatial part of the charge current response is given by 

J\uj) = a'\u)Ej{u) . (15) 

In the language of the previous paragraphs, 0i(x) is an external vector potential A^{x) and the 
operator 0'^{x) is a conserved current J^{x). Choosing a gauge where At = 0, the electric field 



^ This lecture owes a debt to unpublished notes of and private communications with Larry YafFe. 
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becomes = —dtA^. Making a Fourier decomposition where Ax ~ e~^^, we see that = iuA^. 
Comparing eqs. (13) and (15), one sees that the conductivity and the current-current correlation 
function are proportional: 

^«.(^) = . (16) 

Another important player in our discussion will be the heat current Q'^ = T"^ — jjlJ^ , where T*^^ 
is the stress-energy tensor and ji is the chemical potential. We will explain presently why the heat 
current is this peculiar linear combination of the charge and momentum densities. For the moment, 
the important point is that the heat conductivity can be related to two-point correlation functions 
of the stress-energy tensor. The stress tensor couples naturally to fluctuations in the metric g^y. 
Thus, if C'*(x) from above is taken to be T'^^{x), then the corresponding external potential should 
be 5g^u{x). 

Note there is a real and immediate limitation in choosing 0i(x) to be an external vector potential 
or metric fluctuation. By making An{x) external, we are treating the electromagnetic field as a 
control parameter for the system. The system is not allowed to source its own electromagnetic 
fields in this limit, and the photon is not dynamical. Similarly, by choosing 4>i{x) = Sgf^u{x), we are 
forced to work in a fixed background space-time with a non-dynamical graviton. In the limit where 
Coulombic and gravitational interactions are weak, this method of approach makes a lot of sense. 
Indeed, for condensed matter systems, it seems very reasonable to ignore gravity. A non-dynamical 
photon, on the other hand, is questionable; the interactions between the electrons and nuclei are 
electromagnetic in nature. 

Nevertheless, there arc many condensed matter systems in which the photon is treated, to 
first approximation, as non-dynamical. For most metals the Coulombic interaction between the 
electrons is largely screened, and a free fermion gas is a good approximation. In the BCS theory 
of superconductivity, the interaction between the electrons is treated phenomenologically, and 
then the conductivity is calculated from looking at the electric field as an external potential. In 
strongly interacting systems, this electromagnetic nature of the interaction is often obscured and 
renormalized. Our point of view with the holographic models we consider later is a conjecture that 
the strong interactions between the constituents of the system can be separated from the relatively 
weak interaction between the system and an external field strength. 

I would now like to explain why the heat current takes the form Q'^ = T'^^ — ^iJ^ . First, 
consider the case where the chemical potential = 0. The heat current should be the response of 
the system to a temperature gradient. In Euclidean signature, the time component of the metric 
has periodicity 1 /T where T is the temperature. Let the time-time component of the metric have 
the form 

500 = ^ , (17) 

where T{x) is a slowly varying function of position and T(0) = Tq. Euclidean time runs from 
< r < Tq. I assume a possible time dependence of the form diT ~ e"*"^"^. The gradient of yoo is 
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thus the temperature gradient: 



Oigoo - -2^— -2— . (18j 



Consider a change of coordinates in which goo is constant and the temperature gradient is exhib- 
ited instead by a fluctuation of an off-diagonal component of the metric. Under an infinitesimal 
coordinate transformation x'^ ^ x'^ + ^'^j the metric changes by^ 

Sgixu = df^S,^ + d,,ii^ . (19) 

We would like to choose such that di{goo + ^goo) = 0. Setting = 0, we find that 

^9oi = diCo = -^ . (20) 
Continuing back to Lorentzian signature, from eq. (13), we have that 

(TO^-) = GT{u;,0)Sgo.{.) = _ ^r(^.0) ^ . (21) 

Given this result, it makes sense in the absence of chemical potential to identify the heat conduc- 
tivity as 

««M = &^. 

Note that a positive Re[K^'^'] corresponds to a flow from a hot region to a colder region. In- 
deed it must be that Re[K^^] > because of the positivity properties of the spectral density, 
Cc;Im[Gi?(u;,k)] > 0. 

The chemical potential /x is usually deflned as a Lagrange multiplier that introduces an average 
nonzero charge density: H ^ H — jjlQ where Q = J d'^x J* is the charge associated with the 
conserved current J^. Prom this point of view, there is no difference between adding a chemical 
potential /x and introducing a constant background value for the time component of the external 
vector potential At = jj,. 

Given that at nonzero chemical potential, 7^ 0, we have to be more careful with the coor- 
dinate transformation used above. Under such a coordinate transformation, the vector potential 
transforms as well, 

<5^^ = ^,r,^ + ^M,-r, (23) 

which reduces to 6Ai = —AfdiS^o = —jidiT/itoT. 

The change in the Hamiltonian introduced with the external metric fluctuations and vector 
potential is 

5H = - j d'^x {Ti'^dg^, + J^'A^) . (24) 



* The difference between eq. (19) and the Lie derivative (44) is higher order in the temperature gradient. 
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If we are interested specifically in the response of the system to a temperature gradient and an 
electric field, then this change in the Hamiltonian can be written in the form 

^d = -/A((r«-,jO|^w^^) . (25) 

This last equation makes clear that the heat current has the form = T^^ — jiJK 

The linear response of a system to a temperature gradient and an electric field can be summa- 
rized with the following matrix of transport coefficients 

(26) 

Two new conductivities, the thermoelectric coefficients a and a have been introduced above. From 
our discussion above, the conductivities cj, a, a, and k can be directly expressed in terms of the 
appropriate two-point correlation functions of and Q^. In particular 

a=-Gji{u,Q), «= J-GJQ(a;,0) , (27) 

%ix) %Lx) J- 

a=^G^\u:,Q), /c=^GgQ(a;,0). (28) 





B. The Role of Discrete Symmetries 



The transport coefficients a and a are related under time reversal symmetry O. Up to sending 
the magnetic field B —B, the equilibrium state is taken to be invariant under 9. Assume we 
have two real operators ^ and ijj such that 909 = ry^^ and 9^9 = rj^ip where r]^ and rj^ are ±1, 
depending on the choice of operator. The operators T^' and J* are both odd under 9 because 
they are both currents; changing the direction of time changes the spatial direction of the currents. 
Since 9 is an anti-linear operator, we have 

{-BMt,x),i;{0)]\-B) = {B\e[cP{t,x),ij{0)]e\B)* 

= 7i^v^{BM0),(l,{-t,x)]\B) 

where in the last step we used translation invariance. Now if we consider the k^^ = {u>, 0) component 
of the corresponding Fourier transformed, retarded Green's function of cf) and ■0, G'^{u,B), we 
find that 

G^'^(u;,i?) = r?^,?^G^'^a;,-i?) . (29) 
This relation is the Onsager reciprocal relation [18]. The Onsager relation thus implies 



g'^{B) = a^'\-B) , a'^B) = a^\-B) , R'\B) = R^\-B) 



(30) 



15 



In systems with more symmetry, we can put further constraints on these transport coefficients. 
For example, consider a system in two spatial dimensions with rotation and reflection symmetry 
in the presence of a constant external magnetic field. We assume that the equilibrium state of 
our 2+1 dimensional material is invariant under the 90 degree rotation which sends x —y and 
y ^ X. This action of S0(2) demonstrates that the matrices a, a, a, and R, which we collectively 
refer to as M, have the property that M^^ = M^^ and M^^ = — M^^. It is important here that 
we are looking at the k^^ = (w,0) component only of the Fourier transformed Green's functions. 
In the Fourier transform, the integral over cPx washes out the dependence on the rotated spatial 
coordinates. 

Next, we assume our equilibrium state is invariant under a reflection symmetry which sends 
y —y. This reflection sends F^y = B ^ —B. It leaves the diagonal elements of our matrices M 
invariant while multiplying the off-diagonal elements by —1. Together, the 90 degree rotation and 
reflection imply that the matrices cr, a, a, and R all have the property that 

M{B) = M'^{-B) . (31) 

In the case of a and R, this relation is equivalent to the Onsager relation (30). However, for a and 
a, putting this relation together with the Onsager relation, we learn something new, namely that 

a{B) = oF{-B) = a{B) . (32) 

C. The Nernst Effect in High- Tc Superconductors 

Having dwelt long on formalism, it is time to take a break and motivate why these ideas are 

important for the task at hand of understanding condensed matter systems. There are related 
physical quantities which can be derived from eq. (26). For example, there is the canonical heat 
conductivity k, as opposed to R introduced above, which is defined as the heat current response to 
a temperature gradient in the absence of a charge current, 

K = R-Ta-a-^ -a . (33) 

There is also the Nernst effect in which a constant voltage is produced in response to a temperature 
gradient in the absence of a charge current, J* = 0. A little algebra shows that the Nernst effect 
is characterized by the transport matrix 

e = -a-^ ■ a . (34) 

A large Nernst effect has been associated with the onset of superconductivity in high Tg super- 
conductors. For example Wang, Li, and Ong [19] have measured a sharp rise in the Nernst effect 
in the under doped region of La2-a;Srj;Cu04 at temperatures T > Tc. We reproduce a plot from 
their paper as Figure 5. 

In a typical metal, due to an effect called Sondheimer cancellation, the Nernst effect should be 
vanishingly small. Consider a metal at nonzero temperature where the density of states does not 
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FIG. 5: A contour plot for the large Nernst effect measured in La2-:ESr2:Cu04. The Nernst coefficient 
= Oxy/B is plotted in units of nV / K T. Tonset is defined as the temperature at which u begins to differ 
substantially from its high temperature behavior. In the dark blue region, the material is superconducting. 
This figure is described in ref. [19] in more detail. 

vary very much near the Fermi surface. In the presence of a temperature gradient, there are more 
excited electrons above the Fermi surface and holes below the surface in the warmer regions. We 
expect an entropic flux of hot electrons above the Fermi surface to colder regions and a similar 
entropic flux of cold electrons below the Fermi surface to warmer regions. Now if we apply a 
magnetic field, the hot and cold electrons will be deflected in opposite directions. Provided the 
density of states is nearly constant near the Fermi surface, no voltage will be generated, and the 
Nernst effect should vanish. For gold, at room temperature, the Nernst coefficient is on the order 
10 nV / K T [20], compared with values of order 10^ nV / K T measured for La2_xSr2;Cu04. 

Why might under doped, high Tc superconductors have a large Nernst effect associated with 
the superconducting phase transition? Probably because the effective degrees of freedom are not 
electrons and holes. For example, if the effective degrees of freedom are vortices, the Nernst effect 
can be much larger. Magnetic fields penetrate a superconductor through vortex cores, and thus the 
net number of vortices should be proportional to H. Vortex cores have entropy which cause them 
to move in response to a thermal gradient. The phase of the condensate winds around the vortex, 
and thus a vortex moving past a line causes a phase slip of 27r. This phase slippage then induces 
a voltage by gauge invariance. Derivatives of the phase must always appear in combination with 
the vector potential, (9^(/) + A^, and a time varying is gauge equivalent to a voltage. 

The techniques outlined at the beginning of the lecture suggest that we might be able to 
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calculate the Nernst effect without regard to weakly coupled degrees of freedom, be they vortices, 
electrons, or holes. All we need is the two-point functions of the stress tensor and the charge 
current. Remarkably, AdS/CFT provides us with a tool to calculate these correlation functions 
without reference to a weakly coupled microscopic description. More precisely, the weakly coupled 
description provided by AdS/CFT is a dual classical gravitational one in which the field theory 
variables are scrambled almost beyond recognition. AdS/CFT provides a framework for calculating 
in theories where there may not be effective weakly coupled quasiparticles. 

Before turning in the third lecture to a more detailed description of these holographic techniques, 
I first want to get more mileage out of the field theory framework. In particular I want to use 
Lorentz and gauge invariance to place constraints on the form of the current and stress tensor 
two-point functions. These constraints are the Ward identities. 



D. Ward Identities 



In 2+1 dimensions, ref. [21] noticed the following peculiar set of identities relating cj, a, and k 
in the presence of a chemical potential /x and a constant background magnetic field B: 

±a±Tu = (B T /xa;)(7± - (n) , (35) 
±K±Tu = (B T /xo;) a±T - (e + p - /xn) . (36) 

In this expression, I have introduced the energy density T^" = e, the pressure T** = p, and the 
charge density = n. I have also introduced the linear combinations M± = M^y ± iM^x of the 
coefficients of the transport matrices a, a, and R. The remarkable fact about these relations (35) 
and (36) is that if we know a±{io), we know a±{Lo) and R±{lo) as well. There is only one independent 
set of transport coefficients. These relations are really Ward identities, as I now explain. 
Start with a generating functional for Euclidean time ordered correlation functions: 

^wlg,A] = z[g, A] = j D(/)e-^[^'»'^l , (37) 

where the metric g^i, and the vector potential are external and non-dynamical. Given this 
generating functional, we can define the one-point correlation functions as functional derivatives 
of W: 

5W 5W 

Note we have defined T^^"{x) without the customary factor of ^^dci{g^y). Thus T^'^ is a tensor 
density rather than a tensor field. There are similar formulae for the time ordered Euclidean 
two-point functions: 

G^'^'"^(x,y) = {%{T^%x)T^^{y))) = A f ^^^^ , (39) 



aT%,y) . (r.(T-(x)7^(«))) = ^jj^^ . (40) 
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and 

(X. .) . {T.m.mym = J^^^^ ■ (41) 

In the absence of gravitational and gauge anomalies, ^[5,^] should be invariant under diffeo- 
morphisms and gauge transformations: 

x^" ^ x^' + e , (42) 
^ Ai, + d^f . (43) 

Under diffeomorphisms, the change in the metric and vector potential can be expressed as a Lie 
derivative with respect to the vector field 

Sg,.u = {Ci-g)^^ = g^x^^^^ + gu\i^,^ + g,xv,\i^ , (44) 
6A^ = {C^A)^ = + (45) 

That Vrfgi, ^] is diffeomorphism invariant means that 

A short exercise involving integration by parts then shows that 

g,xDf.{T'^''{x)) - Fx^mx)) = , (47) 
where the operator acting on {T^'^) in this expression is the covariant derivative: 

D.iT'^'^ix)) = d.iT^-'ix)) + r;:,(r^''(x)) . (48) 

In fiat space, the Christoffel symbols vanish, and we recover the result 

d^iT^^'^ix)) = F^iJf'ix)) . (49) 

This equation is essentially the Lorentz force law. A much simpler calculation invoking the invari- 
ance of A] under gauge transformations implies current conservation: 

d^iJf'ix)) = . (50) 

The two relations (49) and (50) are the Ward identities for one-point correlation functions. Ad- 
mittedly, we could have derived them more simply from other considerations. However, they are 
the starting point for determining the Ward identities for two-point functions which will in turn 
imply the constraints (35) and (36) on the transport coefficients. 

Two interesting Ward identities for the two-point functions arc obtained by taking a functional 
derivative of eq. (47) with respect to either A^ or Sgfj,i,. I omit the details. The results are 

= ^(%(J"(2/)r^'^(x))) + F/(r,(J°(y) J^(x))) - ^S{x - y)6^''{r{y)) 

+ A-S{x-y)S"'^{J^{y)) , (51) 

= ((r,(r°^(y)T'^'^(x))) + Six - y){g''''T^i'{y) + g^^'T'^'^iy) - g^'^T'^^iy))) 

+Six - y)gP^D^{Tt'%x)) + 5{x - y)g'''' D^{T''P{x)) + F/(T,(r"^(y) J'^(x))) . (52) 
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The covariant derivative is with respect to the x coordinate. I will assume the external field is 

constant in space and time. Assuming translation invariance is not spontaneously broken in the 
equilibrium state, in the flat space limit the one point functions {J^) and {T^^) should be constant 
in space-time. In momentum space, these Ward identities become simpler to write down: 



where I used the one-point function identity (49) to simplify the second relation. We continue 
back to Minkowski space noting that in the scalar case, the analytic continuation of the Fourier 
transformed Euclidean Green's function is the retarded Green's function. For these tensor valued 
Green's functions, we also redefine tensor components with zero indices by multiplying by a suitable 
power of i. Each upper index zero carries a factor of i while each lower index zero carries a factor 



of -i. Thus, we have J% iJ^, Tf iT^, -> -T^, k§ -ik^ = ico, and a similar rule 
for F^'^. With this prescription, the Ward identities become 



From the Fourier transformed Ward identities (55) and (56) and the identifications (27) and 
(28), we can derive the relations (35) and (36). In particular, we consider = (aj,0) in 2+1 
space-time dimensions where we apply a constant external magnetic field F^-y = B and a chemical 
potential At = fi. We assume that the equilibrium state of the theory involves a nonzero charge 
density (J^^) = ((n), 0) and a diagonal stress tensor. 



= -fe;,G^'^'^(fc)-zF/G^''^(A;) + A:'^(J")-M"'^(J'^) , 
= k^ (cf^'^'ik) + S'"'{T'^^') + ^/^^'(T"'^) - S^"'{T''^)^ 



(53) 



(54) 




(55) 



(56) 




(57) 



with constant expectation values for the energy density (e) and pressure {p). As the arguments are 
a little subtle, we present the details in the appendix. 
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III. BEYOND FIELD THEORY: TRANSPORT COEFFICIENTS FROM HOLOGRAPHY 

In this lecture, I will make use of a holographic gravity dual for a strongly interacting 2+1 
dimensional field theory. The field theory in question is maximally SUSY SU(A^) Yang-Mills 
theory in 2+1 dimensions. In 2+1 dimensions, the coupling dimension of mass. If we 

calculate a scattering amplitude or correlation function at some energy scale E, *° appear 

in the dimensionless ratio g\y^/E, and we expect the field theory to be effectively strongly coupled 
at low energy scales. Indeed, there is a widely accepted belief that this theory flows under the 
renormalization group to a nontrivial strongly interacting fixed point at low energy scales. In the 
language of the first lecture, this fixed point is a quantum critical point. Moreover, there is a ten 
year old conjecture — an extension of the original AdS/CFT correspondence — that the infrared 
fixed point is dual to M-theory in an AdS\ x S'^ background [22]. 

We can describe a sector of this conformal field theory with the following four dimensional 
effective action:^ 

S=^J d^x^{R -2A)-^J ct'x^gFABF^'' , (58) 

where A = — S/L^. Consistent with the negative cosmological constant, I will assume that the space 
is asymptotically ^^5*4. Reminiscent of the set-up presented in the second lecture, the Einstein- 
Hilbert term in this action will allow us to compute correlation functions of the stress-energy tensor 
in the field theory while the Maxwell term allows us to compute correlation functions involving 
a conserved current. Conserved currents in field theories come from symmetries via Noether's 
theorem, and the symmetry I have in mind is a U(l) subgroup of the global S0(8) R-symmetry 
present at the infrared fixed point. 
A couple of comments are in order: 

• This classical gravitational description is valid at large where l//t^ ~ N^/"^ . As A^ becomes 
smaller, quantum gravitational effects will become important. 

• Although I seem to have implied that g and k, are independently tunable parameters in the 
action (58), in fact the R-symmetry current is in the same SUSY multiplet as the stress- 
energy tensor, and n and g are related via SUSY [23]: 

= 2g'^L^ . (59) 

• This action (58) describes a sector not only of maximally SUSY SU(A^) Yang-Mills theory 
but of a larger class of strongly interacting 2+1 dimensional SUSY field theories. If I replace 
S"^ with a seven-dimensional object X-j called a Sasaki-Einstein manifold, then there is a 
more general conjecture that M-theory on AdSi^ x X-j is dual to a 2+1 dimensional field 
theory with at least N = 2 SUSY and at least a U(l) R-symmetry (see for example ref. 
[24]). 



^ I will use A, B, C, . . . to index directions in gravity and /i, A, . . . to index directions in field theory. 
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This action S is meant to remind you of the W[g, A] of the previous lecture. RecaUing that 

anti-dc Sitter space has a boundary, the AdS/CFT dictionary says that S, evaluated for a classical 
solution to the gravitational equations of motion, is a generating functional for correlation functions 
in the field theory with the boundary values of g^i, and playing the role of the external metric 
and external gauge field of Lecture II. Thus, S provides a way to compute two-point correlation 
functions of J^" and T^^. 

An important classical solution to the action (58) is a dyonic black hole, i.e. a black hole with 
both electric and magnetic charge: 

^ = 1 {-fiz)dt^ + dx- + dy^) + ~y (60) 
A=—dy-q(l--)dt, f{z) = l + {h^ + q^)a-^ - (1 + {h" + q^)af^ , (61) 

Zh \ Zh/ Zf^ z^ 

where 

(62) 

The radial coordinate z runs from the boundary at z = to the black hole horizon at z = z^. 
Note that the metric approaches that of anti-de Sitter space with radius of curvature L in the 
Poincare patch in the limit z — > 0. In solving for A, a constant of integration was chosen such that 
■^t{zh) = in order that A be well defined at the horizon, AaAbq"^^ < oo. 

Black hole, with its denotation of a spherical object, is a misnomer here. The horizon of this 
black hole is fiat and translationally invariant. Black membrane would be more precise terminology. 
Along with this refinement, it would be more precise to speak of magnetic and electric charge 
density. 

This classical solution is important because it is dual to our strongly interacting field theory at 
nonzero temperature T, magnetic field B, and charge density n. By computing two-point corre- 
lation functions of the stress-tensor and charge current from this background, we are computing, 
by the linear response technology developed in Lecture II, the transport coefficients as a function 
of T, i?, and n. The temperature of the field theory, via the AdS/CFT dictionary, is the Hawking 
temperature of the black hole 



{h? + q-ja 



2^.. 

(63) 



4TTZh 

The magnetic field of the field theory is the boundary limit of the bulk magnetic field Fxy' 

B = h/zh . (64) 

I would like to begin by demonstrating how to use the AdS/CFT formalism to compute (J^) 
given this black hole solution. Near the boundary, the equation of motion for A^ has the solution 



A^ = a^ + h^j,z + . 



(65) 
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On shell, i.e. evaluated for a solution to the classical equations of motion, the Maxwell part of the 
action reduces to a boundary term of the form 



9 J 2=0 9 J 

(We work in the radial gauge ^2 = 0.) Thus the one-point function for the current is 



(66) 



5S 1 



For our dyonic black hole, bt = q/zh and at = —q = fJ,. Thus 

in) = (J*) = = . (68) 

9 Zh 9 Zh 

The electric field of the black hole is thus reinterpreted as a charge density n in the field theory. 
Note that to have a well defined variational problem, the boundary value must be a gauge 
invariant quantity; gauge transformations cannot cause the boundary value to fluctuate. 

As a somewhat more complicated example. Ohm's Law takes an interesting dual holographic 
form. Recall that 

where I am persisting in using these linear combinations of axx and axy from Lecture II and I have 
defined 

J± = Jx± iJy , E± = Ex± lEy . (70) 

We just saw that 

g^J±) = lim d,A± . (71) 

Z—^0 

We can think of ±idzA± as a bulk magnetic field B± and E± as the boundary limit of a bulk 
electric field S±. Prom a gravity point of view. Ohm's Law can thus be written in the form 

a± = lim ^ . (72) 

At this point, I remind the reader about (classical) electric-magnetic duality. For an action of 
the form 

^Jd^xV^FABF^^, (73) 

either Fab or —^eAsCDF^^ could be the fundamental field strength.^ The action is (classically) 
invariant under switching the electric and magnetic fields. For the dyonic black hole background, 
the duality transformation is 

B±^-S±, S±^B±, h^-q, q^h. (74) 



® We define eoi23 
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To compute (7±{n, B) as a function of the charge density and magnetic field requires numerics, 
but there is a constraint from electric-magnetic duality, 

A density plot of |(T+| for various values of the magnetic field and charge density is provided as 
Figure 6. The curious pattern of zeroes and poles bears the imprint of electric-magnetic duality. 
In particular, the duality maps Figure 6a to the negative of Figure 6c and flips Figure 6b about 
the imaginary uj axisJ 

This constraint (75) leads to a very simple way to calculate the conductivity for the field theory 
dual to the uncharged black hole, where h = q = 0. In this case, the background is left invariant 
by the duality, and the conductivity should not change, 

-Mb)?- 



Solving for the conductivity yields 



or in components 



c7± = ±^ (77) 



(^xx ^ g2 ' (^xy = ■ (78) 

This frequency independent result (77) for the conductivity in the absence of a magnetic field 
and charge density, B = p = 0, is surprising. In general, the conductivity could be a function 
of the dimensionless ratio lo/T. This frequency independence was noted originally in ref. [25]. In 
connection with this result, the following quote from ref. [9], which predates ref. [25] by a decade, is 
remarkable: "The distinct physical interpretations of [the high frequency] and [low frequency limits 
of the conductivity] make it clear that, in general, there is no reason for them to have equal values 
(we cannot, of course, rule out the existence of exotic models or symmetries that may cause these 
two to be equal). Apparently, the infrared fixed point of maximally SUSY S'U{N) Yang-Mills 
theory in 2-1-1 dimensions is such an exotic model. 

For the case of nonzero charge density and magnetic field, there is a limit — the hydrodynamic 
limit — in which the conductivity takes a simple, analytic form. The hydrodynamic modes are 
associated with time and distance scales which arc long compared to any microscopic scales but 
short compared to the size of the system. Our system is effectively infinite, and so only the first 
condition is important; given the conformal nature of the underlying field theory, the microscopic 
time and distance scales are set by 1/T. In the limit w/T <C 1 and B/T'^ <C 1, the conductivity 
has the form 

cr± = ±iaQ — (79) 



^ One may also compute the dependence of the current-current correlation functions on k. See ref. [25] for comments 

about the coUisionless to hydrodynamic cross-over in these correlators as a function of k in the case (n) =0. 
® I would like to thank Pavel Kovtun for bringing my attention to this quote. 
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FIG. 6: A density plot of |(7+| as a function of complex u. White areas are large in magnitude and correspond 
to poles while dark areas are zeroes of cr+: a) /i = and q = l,h) h = q = l/V^, c) h = 1 and q = 0. (The 
plots were taken from ref. [26].) 



where 



LOr. 



Bp 
e + p 



7 



e + p 



(80) 



Ref. [26] used AdS/CFT to confirm the result (79). ^ 

Current conservation and the gradient expansion predetermine the form of the two-point cor- 
relation functions in the hydrodynamic limit. More specifically, the facts that d^T'^" = and 
5^ = and the assumption of a well defined gradient expansion for the response of the system 
to small perturbations in the charge, energy, and momentum densities predetermine the form of 
a±, LOc, and 7. The authors of ref. [28] used hydrodynamics to show that the conductivity has to 
take the form (79). 

The gravity dual can then tell us the transport coefficients, for example the normalization of 
the conductivity is 



{e + p^g'^ 



(81) 



Gravity also gives the equation of state, e = 2p, but one could argue this equation was forced on 
us by the conformal scaling which implies the vanishing of the trace of the stress-energy tensor, 
{Tjt) = 0. Away from the hydrodynamic limit, these holographic techniques are arguably more 
powerful, allowing us to calculate for example Figure 6. 



A. Experimental Applications 

Notice the existence of a pole in cq. (79) at the cyclotron frequency. This cyclotron frequency 
LOc is the relativistic hydrodynamic analog of the free particle result cuf = eB /mc. The resonance 
in eq. (79), however, is due to a collective fiuid motion rather than to free particles. There is also 



For an investigation of this system where the constraint S/T^ <C 1 is relaxed, see ref. [27]. 
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a damping 7 that naively one could think of as arising from interactions between the counter- 
circulating currents of positive and negative charged components of the fluid. 

Ref. [28] attempted to estimate the value of this cyclotron resonance for La2-a;Sra;Cu04 in the 
under doped, normal region of the phase diagram: 

This result is about 0.035 times the free electron result. To observe such a resonance in the 
conductivity would be spectacular, even if it says little about the relevance of holographic techniques 
for modeling high superconductors. It says little because the location of the cyclotron pole (80) 
was forced on us by hydrodynamics. On the other hand, observing such a pole suggests that it 
may be reasonable to model the under doped, normal region of La2-xSra;Cu04 using a Lorentz 
invariant field theory where the most important scales are T, B, and n. 

Unfortunately, because real world materials have impurities, this proposal to search for a hydro- 
dynamic cyclotron pole is a difficult enterprise. In addition to providing the estimate (82), ref. [28] 
noted that the inverse scattering time due to impurities in typical samples of La2-a;Srj;Cu04 are of 
the same order of magnitude as the cyclotron frequency, Uc ~ l/Timp- Thus, impurities will likely 
wash out the cyclotron signal. See however ref. [29] for a proposal to search for this resonance in 
clean samples of graphene. 

I demonstrated in Lecture II that a alone determines the transport coefficients a, k, and also 
the related quantities k and 6 using Ward identities. Consider this last matrix, 9 = —a~^ - a, which 
describes the Nernst effect. In the hydrodynamic limit, the Nernst coefficient is 

B ioj 



(83) 



'^xy'-'xx '-^xx'-'xy 

' Xy -•- I •"^cl I J ""c 

Note that in the DC limit (a; — > 0) at nonzero magnetic field {B ^ 0), the Nernst effect 
vanishes, O^y 0. This vanishing is an artifact of translation invariance. As I show in the 
appendix, translation invariance in fact implies that in this limit 

(n) (s) 

^xx ~ ) ^xy ~ ) (^xx ~ ) O^xy ~ • (^^) 

The result for axy is the classical Hall conductivity.^^ 

Clearly, if we arc to describe a large DC Nernst effect, the remedy involves breaking translation 
invariance, and in most materials, dirt and impurities are the most important causes of translational 
symmetry breaking. A proposal was made in ref. [28] to introduce a phenomenological impurity 
scattering time, cj ^ u + i/rimp- With this addition 

Indeed, this result appears to capture some of the qualitative B and T dependence of the Nernst 
effect in high superconductors [28]. 



For the holographic model (58), these results (84) were derived from AdS/CFT in ref. [30]. 
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Sean Hartnoll and I [31] improved slightly on eq. (85) by providing a more explicit holographic 
model of the impurities in which WG were able to capture the 71, and T dependence of Tjjj^p. In 
particular, consider a weak random external potential V{y) coupled to a neutral scalar operator in 
the field theory, 

SH = J dFyV{y)d{t,y) . (86) 
This potential V{y) breaks translational invariance, and we treated it statistically, assuming that 

{V{x)) = , {V{x)V{y)) = VH^^\x - y) . (87) 

The impurities should be a relevant perturbation, important at long distance scales. This condition 
is sometimes known as the Harris criterion: 

[y] = 2 - Ao > , (88) 

where lS.o is the conformal scaling dimension of O. 

Given that V is assumed to be small, we computed the leading order contribution to the impurity 
scattering time using something called the memory function formalism [32] : 

1 y2 /• d^k , 2lmGR'^(a;,k) 
Timp 2x0 '^^0 7 (27r)2 OJ 

In this formula, the momentum susceptibility xo is the long wave-length limit of a momentum- 
density momentum-density correlation function: 

Xo=limG^*'°Xu;,0) = (e + p). (90) 

The result for the impurity scattering time has the scaling form 



But the important point here is given a dual gravity model, we can compute k) exactly. In 

the perturbative limit in V", we have reduced the calculation of Timp to a calculation of a two-point 
correlation function in a gravity model without impurities. (For an attempt to go beyond the weak 
impurity limit using holography and the replica trick, see ref. [33].) 
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IV. HOLOGRAPHIC MODELS OF SUPERCONDUCTIVITY AND SUPERFLUIDITY 

In this lecture, I would like to try to build into the holographic gravity model considered 
in Lecture III the physics of a superconducting or superfluid phase transition. I don't mean a 
quantum phase transition; I mean an ordinary, vanilla, thermal phase transition. We can think of 
the transport coefficients computed in the previous lecture as being the transport coefficients of 
some system in the quantum critical region of the phase diagram, where the physics is determined 
by the effective field theory at the quantum critical point and the most important scales were B, 
/J,, and T. But it would be nice to have, at least in principle, a way of tuning the parameters to go 
through a classical, thermal phase transition — the solid blue lines in Figure 1. 

I propose two methods for modifying the action that will introduce a superconducting or super- 
fluid phase transition. Both methods involve introducing an extra degree of freedom to the gravity 
dual whose boundary value will serve as an order parameter for the phase transition. The first 
proposal is to add a charged scalar field that couples to Fi^^,: 



where D = d — leA [34]. The second proposal is to promote F^^, to an SU(2) gauge field F^^ [35]. 
The order parameter will be a component of the boundary value of A^. 

I would like to make a few general comments about these two proposals before proceeding. 

• It is not, as of the writing of these lectures, totally clear how to embed either of these pro- 
posals in string theory. Hence, we should treat both resulting actions as phenomenological. 
(See however refs. [36, 37] for progress in this direction.) 

• Without a string theory embedding that would fix the potential, y(|^|) is arbitrary and the 
physics may depend sensitively on our choice. We will choose the potential to be a mass 
term, FdV'l) = — 2|'i/'P/-^^^. Note that while the mass = —Ijl? is tachyonic, mildly 
tachyonic scalars are allowed in anti-de Sitter space provided they have masses larger than 
the Breitenlohner-Preedman bound, vr?^p = — 9/4L^ in AdS4. 

• The SU(2) action is fixed by gauge invariance, but the physics is messy. The order parameter 
A'^ is a vector, and the phase transition breaks rotational symmetry. On the other hand, 
one nice thing about the SU(2) case for AdS^ is that the behavior near the phase transition 
is analytic [38, 39]. The symmetries and order parameter of this SU(2) model bear some 
resemblance to those of p-wave superconductors and superfluid helium-3. 

• Normally superfluid phase transitions are associated with spontaneous symmetry breaking 
while superconducting phase transitions with the Higgs mechanism. Thus, we would ap- 
pear to have superfluidity in the boundary field theory and superconductivity in the bulk 
gravity. However, I would like to argue that for some questions, for example in computing 
the conductivity, the difference between a superfluid and superconducting phase transition 
in field theory is not important, and we can pretend the global symmetry group is weakly 




(92) 
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gauged. Understanding the corrections introduced by treating the photon dynamically is an 
interesting direction for future work. 

• The reader should object that there are no classical thermal phase transition in 2+1 di- 
mensional field theories. Infrared fluctuations should destroy any long range order and only 
Kosterlitz-Thouless type transitions are allowed. These models get around the objection 
by having a large number N of degrees of freedom. Given the models' close relationship 
to the maximally supersymmetric S\J{N) Yang-Mills field theory discussed in Lecture III, 
there must be a parameter N that corresponds loosely to the number of colors. I say loosely 
because given the phenomenological nature of the scalar or non-abelian gauge field additions, 
it is no longer entirely clear exactly what N is. 

For simplicity, in the rest of this lecture, I will focus on the scalar case. (More details about 
the SU(2) case can be found in refs. [40, 41].) The dyonic black hole background (60) constitutes 
the normal, high temperature phase. I will return to magnetic fields and the Meissner effect at 
the end of the lecture, but for the moment, I set h = 0. The central observation is that the charge 
density acts as an effective negative contribution to the mass of the scalar [34] : 

mis = m' + g^'Al = m'-'j^,{l- zjznf . (93) 

Note the correction to the mass vanishes at both the horizon and the boundary. However, when g 
is large enough, there will be a region in the interval ^ < z < where m^g becomes too negative. 
There will be a corresponding instability, and the scalar will develop a nontrivial profile. One 
amusing aspect of this model is that there is no need for a term to stabilize the run away 
direction for ^. The curvature of the geometry allows m^jj m? at z = and z = Zh and thus 
stabilizes the run away direction. 

To proceed, I will make a simplifying assumption that <C g^L^- In this weak gravity or probe 
limit, the gauge and scalar sector of the theory do not have enough energy to curve space-time, 
and I can work in a fixed background. Moreover in this limit, the fixed background geometry (60) 
reduces to a black hole metric with f{z) = l — z^/zf^. To further simplify the equations, I will work 
in units where L = e = 1. The equations of motion for the scalar and gauge field become^^ 



z^ 



H'^(f-f)- 

A'l = . (95) 

To have a well defined set of differential equations, I need to specify boundary conditions. The 
physical boundary conditions are as follows. At the boundary {z = 0),ijj and At have the expansion 

At = fi-g'^nz+... , (96) 
= az + bz^ . (97) 



I have made the gauge choice that ip is real. 
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To be consistent with ref . [42] I choose to work in the canonical ensemble where the charge density n 
is fixed. The chemical potential /i is then determined dynamically through the differential equations 
and the other boundary conditions. 

For the scalar, normally the leading behavior would correspond to an external control parameter 
for the field theory, i.e. a source, and the subleading behavior would be an expectation value, just 
as occurred for and the grand canonical ensemble. However, here we have a choice of interpre- 
tations and a corresponding choice of boundary condition [43]. In the AdS/CFT correspondence, a 
bulk scalar with mass m corresponds to a scalar operator in the field theory with conformal dimen- 
sion A through the relation rri^L^ = A(A — 3) (for AdS'4). Thus a value of the mass m^L^ = —2 
corresponds to A = 1 or 2. For a scalar with A = 2, a is interpreted as a source for the operator in 
the field theory while b ~ {O2) is an expectation value. To study phase transitions, we should look 
for solutions with no source for the scalar, a = 0, but where b becomes nonzero at some critical 
temperature. On the other hand, for a scalar with A = 1, things are switched, a ~ {Oi), and one 
should set 6 = 0. With a single choice of mass parameter, I get two models for the price of one. 

The other boundary conditions for this field theory I set at the horizon of the black hole, z = z^. 
For the scalar, the condition -!/; < 00 at the horizon eliminates one of the integration constants in 
the differential equation. On the other hand, for the gauge field, I need to choose Af = at the 
horizon in order to ensure the g^^Af < 00. 

Given the boundary conditions and the set of differential equations for At and tp, I have a well 
posed problem which unfortunately does not appear to have an analytic solution. However, the 
differential equations are relatively straightforward to solve numerically, and one can look for a 
phase transition as a function of the dimensionless ratio n/T^. In practice we fix the temperature 
by setting the horizon radius to = 1, and tune n. For < n < ric, the black hole solution with 
ip = appears to be stable. However, for n > ric, there is a phase transition to a black hole with 
scalar hair. We can equally well think of this phase transition from an ordinary black hole to a 
hairy black hole as occurring as we lower the temperature at fixed n. Plots of the behavior of {Oi) 
as a function of temperature are given in Figure 7. The phase transition is second order (see for 




FIG. 7: The condensate as a function of temperature for the two operators Oi and 02- The condensate 
goes to zero at T = Tc oc n^^^. 
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example the free energy plot in ref. [44]). For T <Tc, the order parameters scale as 

(Oi)~(02)~(re-r)V2 (98) 

which has the classic 1/2 mean field exponent of Landau- Ginzburg theory. 

We will see after studying the conductivity of the system that there is a sense in which wc can 
interpret (Oi) and {O2) as twice the superconducting gap. Recall that in BCS theory, there is a 
classic prediction for this number at T = 0, namely 2 x gap = 3.54 Tc [45]. In comparison, in these 
holographic systems, wc find that 2 x gap is cither infinite in the A = 1 case or about 8Tc in the 
A = 2 case. The relatively large size of the gap compared to BCS theory suggests our holographic 
field theory is strongly interacting. Note that for the high Tc superconductors, 2 x gap is typically 
in the range 4 to 7 times Tc [45].^^ 



A. Conductivity 

As we did before in the dyonic black hole case, we can study the conductivity of the field theory 
dual to this hairy black hole. To calculate the conductivity, we need to solve the equation of motion 
for a fluctuation of the gauge field in this background: 

ifA'x)' - = %i.'Ax . (99) 

The near boundary solution for can be written 

A^ = ^ + ^-rz + ... (100) 

loj a 

where we have assumed that Ex ~ e"*"^*. The longitudinal piece of the conductivity matrix is then 



X 



(Txx=ir ■ (101) 

i^x 

In general, wc again need numerics to solve this equation, but there are two nice limits in which wc 
can study the behavior analytically. For T > Tc, the scalar vanishes tjj = 0, and axx = l/ff^, as wc 
calculated in Lecture III. For T ~ and the scalar with conformal dimension A = 1, numerically 
we find that ip is a nearly linear function of z: 

Given this numerical observation, eq. (99) reduces to the Klein-Gordon equation: 

A'^-u;^Ax = {0^fAx . (103) 



The interpretation of {Oi} in terms of an energy gap is to some extent cheating. The strong interactions imply 
we don't necessarily have a good particle interpretation. Even if we did have a good microscopic picture of a 
condensate of some kind of bound state, we could have n-particle bound states in place of Cooper pairs. 
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Note that at T = 0, the warp factor / = 1 and the geometry extends from z = to ^; = oo. This 
Klein-Gordon equation has two solutions, depending on the magnitude of lo: 

In the first case, we chose the z ^ oo boundary condition by demanding the solution be finite. In 
the second case, we chose the large z boundary condition by demanding the solution be a wave 
traveling to larger z. Thus, the conductivity takes the form 



a.. = \^^^^^s^ {{0,f - o;^) . (105) 

In principle, we should have been more careful about an ie prescription for the w ^ limit 
of eq. (105), but wc can recover the full behavior from considerations of causality and analyticity. 
There must be a delta function peak at u; = by the Kramers-Kronig relations: 

Imcr(a;) = V\ ^ . (106) 

TT LO' -U 

Thus, 

Rec7(a;) ~ 7r(Oi)(5(a;)/52 (107) 

near a; = 0. This formula for a^x produces a graph that looks very much like classic textbook plots 
of the conductivity calculated for BCS superconductors. Note the size of the superconducting gap 
would be (Ci)/2. The factor of one half comes from the fact that in BCS theory the dissipative 
conductive response involves producing a pair of excitations. A plot of this T — limit of the 
conductivity is given in Figure 8. Ref. [42] examined the behavior of the conductivity at interme- 
diate temperatures, < T < Tc. As expected as one increases T, the gap in the real part of the 
conductivity begins to fill in.^^ 

Ref. [48] studied this gravitational Abclian-Higgs model away from the weak gravity limit. I do 
not have the time to describe those results here in detail. There are many qualitative similarities. 
There are also two interesting differences. The first concerns translation invariance. Away from 
the weak gravity limit, the normal component of the fluid is charged, and thus the conductivity a 
is a nontrivial function of w in the normal phase. The second is that there is an additional source 
of instability associated to the fact that the effective Breitenlohner-Preedman bound changes as 
the charge on the black hole increases. Given eq. (93), one might expect that neutral scalars 
will not produce an instability. In fact, as the charge on the black hole approaches extremality, 
the geometry near the horizon approaches AdS2- The BF bound associated with this AdS2 is 
> —3/2, implying that neutral scalars with a mass —9/4 < < —3/2 will condense to 
form hairy black holes at low temperature. From the field theory point of view, this instability is 
associated with a renormalization of the effective potential for the order parameter at large density. 



A Mathematica notebook [46] that reproduces the numerical work described in ref. [42] is available at the website 
[47]. 
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FIG. 8: A plot of the conductivity in the hmit T for the scalar Oi. The solid black line is the Kea^ 
The dashed blue line is the Im axx rescaled by a factor of five. The dotted red line is the conductivity a, 
(in the probe limit) when T > Tc- 



B. The London Equation 



Given the importance of magnetic fields in the phenomenology of superconductors, I would like 
to explain how the London equation arises in this model in the weak gravity limit. The London 
equation, 

J(u;,k) = -nsA{u,k) , (108) 

was proposed (in a gauge where the order parameter is real) to explain both the infinite DC 
conductivity and the Meissner effect of superconductors. This equation is valid in the limit where 
both UJ and k are small compared to the scale at which the system loses its superconductivity. In 
our model, that scale is set by {Oi). One important and subtle issue in understanding this equation 
is that the two limits w — > and k ^ do not always commute. In the limit k = and w — > 0, we 
can take a time derivative of both sides to find 

J(cj,0) = — E(u;,0) (109) 
u> 

explaining the infinite DC conductivity observed in superconductors. On the other hand, in the 
limit UJ = and k — > 0, we can instead consider the curl of the London equation, yielding 

V X J(x) = -n^B(x) . (110) 

Together with Maxwell's equation V x B = 47rJ, this other limit of the London equation implies 
that magnetic field lines are excluded from superconductors: 

-V^B = V X (V X B) = 47rV x J = -47rn,B , (111) 

which has exponentially damped solutions. The magnetic penetration depth squared = I/AttHs 
is proportional to the inverse of the superfluid density n^. 
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An important point is that because the gauge field is external, there is no Maxwell equation 
in the holographic model: Currents in the material to not source electromagnetic fields. Thus we 
only get half of the Meissner effect; we only have the London equation. 

In the discussion of cr, I considered only the first limit, having set k = 0. However, the London 
equation holds more generally, including in the limit where uj is sent to zero first. Including the 
k dependence in eq. (99), the solution (104) can be generalized by replacing u;^ with the Lorentz 
invariant combination — k^. Thus, in the present case the limits a; — >^ and k — > commute 
and one finds 

3{uj, k) «i -^^A(a;, k) (112) 

for small oj and k, allowing us to identify {Oi) / with the superfluid density n^. A similar equation 
holds for the model with A = 2, but the confirmation requires numerics [48]. 

Many groups have studied this holographic model of a superconductor during the past few 
months. Given my restriction above to a 2+1 dimensional field theory and an order parameter 
with A = 1 or 2, I would like to single out ref. [49] for special notice. The authors considered a 
small selection of related models in both 2+1 and 3+1 space-time dimensions with various values 
for A. Their results are qualitatively similar to what we have presented above. Ref. [50] lists more 
related work. 



C. Second Sound 

Despite its more rigorous interpretation as a superfluid, thus far I have discussed this holographic 
model as a superconductor. In this last part of Lecture IV, I would like to address the more standard 
interpretation of the model and discuss a classic phenomenon associated with a superfluid, second 
sound. Second sound is a collective motion available to fluids with two components where the 
components move out of phase with respect to each other. 

As discussed in ref. [44], a superfluid has a pressure that depends on three thermodynamic 
variables, 

P{T,^i,diip) , (113) 

where tp is the phase of of the condensate, = j'^je*'^. The quantity diip is sometimes called a 
superfluid velocity. Turning on diip sources a charge current in the system. Coupling the order 
parameter to an external gauge field, the pressure can be written 

P{T,-Do^,Diip) , (114) 

where D^tp = dfj,ip — A^. As we did before, we are free to make a gauge choice where the phase of 
the condensate vanishes, leaving P a function of A^, 



p(r,Ao,-^i) 



(115) 
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T T 

a) M h) M 

FIG. 9: The speed of second sound as a function of r/|/i| , computed by evaluating thermodynamic derivatives 
in Eq. (116): a) Oi scalar, b) O2 scalar. The speed of second sound vanishes as T Tc and appears to 
approach a constant value as T — s- 0. 



The second sound mode can be derived from a hydrodynamic analysis of the current conservation 
conditions. We have d^T^" = and d^J^ = 0. Working with Afj_ instead of because of the 
equality of mixed partials, we need to impose that doAi = diA^. 

For simplicity, we will continue to work in the weak gravity limit. In this limit, the conservation 
condition for T^^" is irrelevant for the second sound mode. Ref. [44] demonstrated in this case and 
in the absence of a charge current, the second sound speed has the form 



2 _ _ d'^P/d'^Ai 
^2 ~ " ~[PpJ[PAq 



(116) 

A,=0 



As the normal component of the fluid is stationary in this probe limit, ref. [51] argued that there 
is a close analogy between this quantity and what is called fourth sound in the superfluid helium 
literature. In fourth sound, superfluid helium is forced to travel through a capillary that has been 
packed with a powder that only lets the superfluid component move freely. The fourth sound is a 
collective motion of the fluid in such a capillary. 

A plot of the second sound speed versus temperature for our two models with A = 1 and 2 is 
shown in Figure 9. Note that w| ~ (Tc — T) near the phase transition. 
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V. A COMMENT 

I have emphasized high temperature superconductivity in these lectures because it is one of 
the most interesting and important among the unanswered questions in condensed matter physics. 
However, it may prove that high Tg superconductivity is not a good target for these holographic 
techniques. One should keep an open mind. Heavy fermions, fermions at unitarity, transitions 
between fractional quantum Hall states, or yet some other condensed matter system may be more 
suited to an AdS/CFT approach. A dream is to one day find a material whose essential physics is 
well described by a gravity dual and, more than that, where the dual predicts some new physical 
effect. In the meantime, it is instructive to learn about the properties of a class of strongly 
interacting field theories and fascinating to see how transport coefficients and phase transitions in 
field theory can be mapped to black hole physics. 
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APPENDIX A: DERIVING THE WARD IDENTITY CONSTRAINTS ON THE 

CONDUCTIVITIES 

We begin by specializing the Ward identities (55) and (56) to the 2+1 dimensional case of 
interest. For spatially homogeneous fluctuations where k'^ = {lo,0), eq. (55) becomes 

ojG^\ij) - iBekj&^{u) + 5iju{n) = , (Al) 
where we let the indices i, j, and k run over the values x and y. In components we thus have 

u&^'\uj) + iB&j^y{uj) + tj(n) = , 

u&^'\u) - iBGYii^) = . (A2) 

From the Ward identity (56) relating the stress-tensor stress-tensor correlation function to the 
stress-tensor current two point function, we find that 

u&^^\u) + 6ij {e)uj - iBekj&^^{uj) = (A3) 
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or in components that 

a;G^"'°"(L^) + oj{€) + iB^^'^ioj) = , 



tc;(5°"'°^(6c;) - iB&^'^'iLo) = . (A4) 



(JJ)± 






{JT)± 






{TT)± 







At this point, in analogy with the notation used for the transport coefficients in eqs. (35) and 
(36), we find it convenient to introduce a complexified notation where 

(A5) 
(A6) 
(A7) 

With this notation, we can replace our 2x2 antisymmetric matrices of transport coefficients with 
complex numbers. The discussion after eq. (30) implies that {JT)± = {TJ)±. In terms of the 
complexified combinations, we find that eq. (A2) becomes 

±uj{TJ)± - B{JJ)± + co{n) = , (A8) 

while eq. (A4) reduces to 

±u{TT)± - B{TJ)± + (e)a; = . (A9) 
Prom eqs. (27) and (28) it follows that {JJ)± is related in a simple way to a±: 

{JJ)± = iG'^jfiiu) - = {a^y ± ia^^)cj = ua± . (AlO) 

Similarly, one finds that 

{QJ)± = ua±T ; (QQ)± = loR±T . (All) 

We can calculate the u ^ limit of {JJ)± and {JT)± from a Lorentz boost argument. Consider 
a plane with charge density (n) carpeted by a perpendicular magnetic field B. If we boost the 
plane by a velocity v in the x direction, then to linear order in v, we find an electric field Ey = vB 
in the y direction. Moreover, there is a current Jx = {n)v. Thus we find that 

lim axy{ijj) = lim -{JJ)± = {n)/B . (A12) 
o)— >0 w— >0 (jj 

To calculate {TJ), we start with a diagonal stress tensor {T^^) = (e) and (T^^) = {p) in the rest 
frame of the plane. After the boost, we find (T"*) = (e + p)v. Thus 

lim -{TJ)± = {e + p)/B . (A13) 

There is a potentially subtle problem with the w ^ limit of Eqs. (AS) and (A9). A retarded 
Green's function involving a momentum density and a translationally invariant equilibrium state, 
such as {TT)± or {JT)±, should vanish in this limit by translation invariance. However, the Green's 
functions as defined through the Ward identities may have contact term corrections which prevent 
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this vanishing. Plugging eqs. (A12) and (A13) into eq. (A8), we see that {JJ)± and {TJ)± vanish 
in a way that is consistent with the Ward identity. However, there is a problem with eq. (A9): 

lim(rr)± = ^(p) . (A14) 

UJ—^0 

Thus, we should correct eq. (A9) by adding a contact term {p)uj, 

±u;{TT)±- B{TJ)± + u{e + p) = . (A15) 

We can re-express eqs. (A8) and (A15) in terms of instead of r°', yielding 

±oj{JQ)± + {±fiuj - B){JJ)± + Lo{n) = 0, (A16) 
±u{QQ)± + {±i^oj - B) {JQ)± + uj{e + p- nn) = {) . (A17) 

Using the eqs. (AlO) and (All), these eqs. (A16) and (A17) reduce to the relations (35) and (36) 
as promised. In employing these relations, it is helpful to keep in mind that the limits a; — and 
5 — > often do not commute. 
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